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TECHNICAL PAPER 


BUCKLING OF COMPOSITE BEAMS 
INTRODUCTION 


Torsion/bending coupling is an interesting, inherent phenomenon occurring in open-section 
structural elements subjected to bending and/or combined axial/bending loads. The benefits of using 
open-section beams are many: stackability, weight savings, and dual functions (e.g., serving as load 
carriers as well as conduits). Likewise, composites offer advantages over their metallic counter- 
parts; primarily, their higher strength-to-weight ratios and lower coefficients of thermal expansion 
(CTE) are favorable. This type beam is used extensively in the aircraft industry in the fuselage sup- 
port structure and wing supports, and has present and potential future use in aerospace applications; 
namely tankage stiffeners, intertank stiffeners, and as primary load carriers in truss work. However, 
in each of the aforementioned applications the layups are symmetric in nature. By understanding the 
bending/torsion effects indigenous to open-section beams one can make use of the tailorability of 
composites (to achieve desired optimum material properties) without compromise to the design 
(e.g., no added plies to accomplish symmetry). 

Antisymmetric layups lend themselves well to thin-walled structural members. With an 
average ply thickness of 0.005 in, a beam of wall thickness less than 0.040 in (8 plies) is difficult to 
manufacture using a symmetric laminate layup. 

Even an open-section beam made of an isotropic material has an inherent tendency to bend 
and twist if the transverse load is applied anywhere other than along the plane of the shear center. 
A transverse load applied to such open sections undergoes pure bending only when the load is 
applied at points along the shear center located outside the section. 1 As Valsov illustrates, an open 
section not loaded at points along its shear center will twist and bend as in figure 1(b). Figure 1(a) 
shows the same open-section beam with no load applied. In the case with loading at the shear 
center, notice that pure bending results, as in figure 1(c). 

As in symmetrically laminated structural elements, unsymmetrically laminated elements can 
have many different resultant properties depending on the chosen material system, ply orientations, 
and layup sequence. 2 Add these factors to the bending/twisting coupling inherent to an open section 
loaded transversely and one can easily surmise the difficulties that might be associated with design- 
ing using unsymmetric, open-section structural elements. 

Characterization of the relationship between angle of twist and ply orientation is one of the 
objectives of this paper. This is being done based strictly on test results, with results of the analyti- 
cal aspect to be presented in a later paper. Also examined and presented here are the results of 
testing involving shear center location, load application relative to shear center, and the resulting 
angle of twist on open-section beams. Curves depicting these parameters and the various inter- 
relationships are contained in this paper. 




APPROACH 


To reach the first objective (characterization of the relationship between the angle of twist 
and ply orientation), the effects of distortions due to CTE must be included. This produces an amount 
of pretwist in the beam apart from the twisting that is to occur due to the load not being applied at 
the actual shear center. Examination of the various laminate layups will reveal a general relationship 
between the ply orientation and angle of twist. 

Open-section beams of various ply orientations were manufactured and tested to determine 
what effect ply orientation and laminate layup have on the shear center location, i.e., can the shear 
center in open-section structural elements be shifted toward its geometric center? As illustrated in 
figure 2, the geometric center and theoretical shear center differ. Yet this moving of the shear center 
can possibly be observed by comparing the various twist and deflections in each beam tested as a 
specified load is applied at designated positions along the horizontal plane of the shear center 
(fig. 2). Recalling that no twist (pure bending only) will occur when transverse loads are applied at 
the shear center, one can better define the location of the “true” shear center. This “true” position, 
since it is observed after the curing process (hence, CTE mismatches are already accounted for), 
reflects the dependency of shear center on laminate layup. No attempts were made to quantify this 
preload twist; however, measurements were taken prior to testing the beams to analyze trends. 

Open-section beams (C-channels) of the layup [O/0] m are referred to as Form 1. The 
laminates referred to as Form 2 are of ([O/0] n ,[O/0] n ) construction. 



Figure 2. Load application points. 
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MANUFACTURING METHODS AND TESTING 


As previously stated, the beams were fabricated using the hand layup and hot-drape forming 
techniques. Hand layup involved cutting 3-in wide strips of the unidirectional tape to the correct 
lengths and wrapping the tape at proper angles across the tool (in this case a solid aluminum square 
tube served as the tool) layer by layer. Intermediate vacuum bagging was performed to ensure ade- 
quate compaction prior to the autoclave cure of the part. Hot-drape forming was a new procedure 
used to take advantage of the multiaxis tape laying machine. Using this technique panels 
(approximately 18-in wide by 72-in long) were tape layed by machine. This provided repeatability of 
the part, provided adequate compaction (tape head is pressure sensitive, i.e., constant pressure as 
tape is layed in place), and reduced time of manufacturing. After the panels were cut to a strip 9 by 
72 in, they were “draped” over the same tool used previously, bagged, and placed in the hot-drape 
forming “oven” under controlled temperatures (see the photographs in the appendix). This allowed 
the material’s resin to flow enough so that the panel would take the shape of the tool. 

After the part was formed (regardless of method), it was vacuum bagged for autoclave cure 
using standard procedures. The resulting channels were then removed from the tool and cut and 
machined to dimensions suitable for testing. Figure 3 shows the geometry of the finished beam prior 
to instrumentation for testing. 



C-CHANNEL CROSS-SECTION C-CHANNEL COMPOSITE BEAM 


Figure 3. C-channel beams. 



The testing included applying load at several locations along the horizontal plane of the shear 
and geometric centers. After the beam was instrumented and the load cells calibrated, load was 
applied in 15 to 20 lb increments until a total of 130 lb was being applied. Beginning at the theoretical 
shear center, this procedure was repeated (at 0.25-in increments) until the final measurements at 
the geometric center (fig. 2). During the test procedure, the deflections and strains were recorded 
with the addition of load at each of the incremental points. 


EFFECT OF FIBER ORIENTATION OF FLANGE DEFLECTION 


Graphs 1 through 6 of section I of the appendix show deflections in the upper and lower 
flanges for samples of various layups of the C-channels. Also, charts 1 through 6 show a complete 
set of deflection data for the [(0/75)] 12 case. The beams were instrumented as shown in figure 4. An 
interesting occurrence in the beams, discovered through observation of the test results, is that the 
deflections decreased as the load application point is moved away from the center of gravity. This 
fact becomes more important in the following section when angle of twist is discussed. As expected, 
the deflection represented by D9 on the charts is smallest; this measurement is taken at a point on 
the beam’s web. The top flange deflected the most in all cases examined. 



Figure 4. Deflection gauge locations for composite beam load tests. 


Another observation was that for beams having the same ply angle, the beams of Form 1 
deflected at least twice (100 percent) as much as those of Form 2 when the load was applied near 
the shear center, and almost 80 percent more when loaded at the geometric center. Additionally, by 
examining the deflection data at the 0.25-in increments, more interesting facts were revealed. For 
example, in the case where the cross ply angle is ±75°, it is noted that the deflection in the web and 
flanges decrease (for the maximum load case) as one moves toward the calculated shear center 
located 1.44 in from the geometric center. Yet, as will be discussed later, the deflections approach 
zero in the web even more by moving further than 1.44 in from the shear center. This serves as an 
indicator that the true shear center is further than 1.44 in away from the centroid. Therefore, one 
would opt for the laminate layup of Form 2 in situations in which low deflections are desired. 

Two materials systems were used: AS4/3501-6 and IM7/8552. Figures 5 through 8 show 
deflections measured at the geometric centers and shear centers for a [0/75] 12 beam of each 
material. Due to material properties, the IM7/8552 beams deflected more than the AS4/3501-6 
beams. Another point of interest from the data gathered is that the further away from the geometric 
center that the load is applied, the greater the difference in the deflection amounts. The explanation 
of this could also form the basis for more research work. 
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DEFLECTION 





DEFLECTION 






DEFLECTION 


IM7/8552 [0/75] 1 2 



LOAD 

Figure 7. Deflections for C-channel loaded at geometric center (IM7/8552). 


IM7/8552 [0/75]1 2 



LOAD 

Figure 8. Deflections for C-channel loaded at theoretical shear center (IM7/8552). 
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EFFECT OF FIBER ORIENTATION OF SHEAR CENTER LOCATION 


For this research effort, the theoretical shear center for the given cross section is easily 
determined by the formula: 


e = bht/41 . 

After this shear center is known, one would expect that a transverse load applied through 
this point would produce only bending, and that the same load applied through the beam’s geometric 
center produces both bending and torsion coupling. Yet, due to the layup of these open-section 
members (no midplane symmetry) the true shear center does not follow this equation. The research 
work showed that for some of the laminate layups the true shear center actually lies between the 
geometric center and theoretical shear center, and could, in some cases, lie beyond the calculated 
shear center. As the angle increased, the amount of twist observed and measured when the beam 
was loaded at this “predicted” shear center increased also. As the load application point was 
shifted toward the center of gravity (c.g.), the trend was for more twisting to occur. However, for 
most of the cases tested, this increase did not occur until after the load application point had been 
shifted more than 0.25 in from the shear center in the direction of the c.g. This would indicate that the 
true shear center is actually between the predicted shear center and a 0.25 in toward the c.g. Yet, for 
those beams of Form 1 with a cross-ply angle greater than 45°, the “true” shear center was outside 
the theoretical value, i.e., more than 1.44 in from the geometric center. 

Figure 9 shows the angle of twist for various layups as calculated by the formulas 3 

e = d\/2L and e = d 2 /2L , 

0 = di/(h/2) and <p = d 2 /(b/ 2) , 


where 


h = height of web 
b = width of flange 
d\ = web deflection 
d 2 = flange deflection. 

For presentation purposes, beams made with only the 15°, 45°, and 75° cross-ply angles of 
Form 1 are shown. Measured strains and deflections for the 30° and 60° cases fell between the 
extremes of the charts shown. Charts 1 through 6 in section II of the appendix contain the strain data 
as measured for the [(0/1 5)] 12 beam. 

Graphically, one can see from figure 10 that the angle of twist does increase with an increase 
in the cross-ply laminate angle, regardless of the overall layup (Form 1 or Form 2). The range in 
angle of twist as measured is between -18° and 55° for Form 1, while its range was -25° to 75° for 
Form 2 beams. 
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Figure 9. Theoretical versus measured shear center comparison. 


What these equations and formulas do not show, however, is the postcure inward canting of 
the beams as shown in figure 7. No quantitative assessment of the amount of this cant is presented 
here. Nevertheless, it should be pointed out that this phenomena occurred in all the beams. Interest- 
ingly, the beams containing 45° plies of Form 1 canted outward (i.e., pulled away from the tool after 
curing) and twisted along the length of the beam by more than 90°. CTE’s, as well as interlaminar 
shear stress and strains during the curing process, would have to be closely examined. The explana- 
tion of this occurrence could provide the basis for additional research work. 


Test data were also used to gain insight on the location of the beam’s “true shear center 
(i.e., no twist when transverse load applied) as a function of laminate layup. Due to the beams 
geometry, they are expected to twist as well as bend, and all the beams tested experienced this. 
However,’ the focus here was to attempt to categorize the amount of twist that could be expected 
from a given layup. In analyzing the test results, several interesting trends were noticed and are 
reported here. 

First, the beams of the Form 1 construction had higher degrees of rotation (angles of twist) 
than those of Form 2 for any given cross-ply lamina. Remembering that the Form 1 beams have 
cross plies only in the positive direction, whereas the Form 2 beams have both positive and negative 
plies, it is obvious that, upon completion of the appropriate stiffness matrices, this in itself is not 
abnormal. 


Figure 10 reveals another trend. By calculations using the given geometry of the beams, one 
can find that the theoretical shear center is located 1.44 in to the left of the geometric center as 
shown in figure 3. By observing the charts, one can see that none of the plotted beams show an 
experimental shear center that matches its theoretical value. All the beams were of the same overall 
geometry and the results shown on the charts were based on a 130-lb transverse load placed on the 
beams. 
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Figure 10. Angle of twist 


[0/15] Laminate Layup 
Angle of Twist 



[0/45] Laminate Layup 
Angle of Twist 



[0/75] Laminate Layup 
Angle of Twist 



Distance from Geometric 
Center to Shear Center 

sus load application points for Form 1 beam (AS4/3501-6). 


The [0/75] 12 layup most closely approached the theoretical value and had the highest angle of 
twist, while the [0/15h 2 layup had the lowest angle of twist. This tendency, particularly the least 
twist in the [0/15] !2 members, was not expected. It would appear intuitively that since the tendency 
to rotate in these beams would not be offset by the almost horizontal plies of the 15° laminates they 
would have the higher twist angle. Similarly, these same beams had true shear centers located more 
than 1 inch from the theoretical value. 

All beams had some small amount of twist before loading due to CTE differences between the 
plies during the cure process. This pretwist was most noticeable in the [0/45] i 2 beams which rotated 
more than 90° along the length of the beam. Pretwist due to CTE mismatches is not addressed here, 
but it is worthy of future consideration, as is controlling the various parameters comprising the cure 
cycle (e.g., cure pressure, tooling materials, cure temperature profiles). This pretwist, however, does 
not significantly alter the trends reported here since the loads were applied at the calculated theo- 
retical shear center and additional twisting occurred, again verifying the shifting of the “true” shear 
center. Figures 11(a) and (b) shows the angle of twist for Form 1 beams when 0 is equal to 15, 45, 
and 75, respectively, and only a minimal load is applied (P < 5 lb). 

[0/15] Laminate Layup 
Angle of Twist 



Distance from Geometric 
Center to Shear Center 

Figure 11(a). Angle of twist versus load application points for Form 1 beam 
with minimal load applied, P < 5.0 lb (AS4/3501-6). 
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Distance from Geometric 
Center to Shear Center 

Figure 1 1(b). Angle of twist versus load application points for Form 1 beam 
with minimal load applied, P < 5.0 lb (AS4/3501-6). 






Assuming plate theory, it can be shown that in-plane stress resultants for laminated struc- 
tures are not only functions of midplane strains but can be functions of curvatures and twists. Addi- 
tionally, the in-plahe forces can cause deformations that can cause curvatures or twisting deforma- 
tions. 4 Plate theory does not directly apply to beams as presented in the above cases. As Vinson 
points out, the above analytical procedures are valid only if no coupling exists, i.e., the £>i6 and £>26 
terms of the stiffness matrix are zero. 5 The equations were used in this paper only to produce a rela- 
tive number that would allow numerical comparisons to be made. More complex analysis that con- 
siders not only the nonzero cases for £>16 and £>26 but also hygrothermal effects would need to be 
developed to more accurately assess the true values for the amount of twist produced due to lami- 
nate angle orientation. 

Even though plate theory may not be directly applicable to beam theory, one can use it to 
gather trend data. This is an attempt to categorize the behavior of the laminates rather than assign a 
quantitative value for the angle of twist. Regardless of the number assigned to this angle, the 
observed trends would still be valid. The analysis method used to arrive at this value does not influ- 
ence the direction or importance of the trend. 


CONCLUSIONS 


The magnitude of information gathered through tests performed on open-section unsymmetric 
beams was tremendous. All data gathered were not discussed in this paper, however, subsequent 
writings will explore other areas worthy of further investigation in an attempt to more accurately 
categorize the behavior of such beams. Major findings reported here deal with the primary objective 
regarding shear center location and manipulation, as well as angle of twist and warpage in these 
structural members. 

One important revelation is that the angles of twist and associated warpage in open-section 
beams is higher for antisymmetric beams of the [O/0]m laminate layup. Some amount of bend- 
ing/stretching is expected in any antisymmetric layup due to the fact that the components of the [B] 
coupling matrix, which relates stress couples to midsurface planes and stress resultants to curva- 
tures, are not all zero if the structure is not symmetric about its midplane. Likewise, unless the £>16 
and D26 terms of the stiffness matrix [£>], which relates the stress couples to curvatures, are zero, 
bending/twisting will occur. ^ No attempt is made to measure the magnitude of hygrothermal effects 
on twist in the beam, even though there could be as many as three different CTE’s depending on 
cross-ply angle and orientation, all which play a role in the warping of the beams. ^ 

The major fact borne out here is that laminate layup does control shear center location. As 
evidenced by charts presented earlier, the same geometric shape can yield “true” shear centers that 
differ by more than an inch depending on laminate layup. The location of the “true” shear center for 
the type beam investigated here depends on several factors. Cross-ply angle sequence is important 
(Form 1 or Form 2) in that other unsymmetrical layups using the same angle would yield different 
shear center locations. Yet not to be overlooked are factors such as cure cycle and laminate materi- 
als. 
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" 1T1009“ j 

" ITIOIO"^ 

" 1T1011" 

" 1T1012” 

" 1T1013" j 

■ 2T1009" 

* MST " 
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“ MST “ ! 
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-1.04522 

-1.04236 

-1.04527 
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-25.60925 

155.94434 

-8.84251 

-34.41862 
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Manufacturing and Testing Photographs 
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Hot-drape forming "oven." 
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Laminate placed on tool ready for hot-drape forming. 
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SLACK AND WHITE PHOTOGRAPH 


C-channel upon removal from tool after autoclave 
Notice warpage of beam along its length. 
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Fully instrumented beam C-channel prior to testing. 
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Front view of test configuration. Slide block allows load 
application point to be changed. 







Test being performed on L-angle. Beam has started to rotate. 
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